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Abstract 

Bound (/(/-systems are studied in the framework of different 3-dimensional 
relativistic equations derived from the Bethe-Salpeter equation with the in- 
stantaneous kernel in the momentum space. Except the Salpeter equation, 
all these equations have a correct one-body limit when one of the constituent 
quark masses tends to infinity. The spin structure of the confining qq inter- 
action potential is taken in the form 2:7172 + (1 — x)I\l2, with < x < 1. 
At first stage, the one-gluon-exchange potential is neglected and the confin- 
ing potential is taken in the oscillator form. For the systems (us), (cu), (cs) 
and (uu), (ss) a comparative qualitative analysis of these equations is carried 
out for different values of the mixing parameter x and the confining potential 
strength parameter. We investigate: 1) the existence/nonexistence of sta- 
ble solutions of these equations; 2) the parameter dependence of the general 
structure of the meson mass spectum and leptonic decay constants of pseu- 
doscalar and vector mesons. It is demonstrated that none of the 3-dimensional 
equations considered in the present paper does simultaneously describe even 
general qualitative features of the whole mass spectrum of qq systems. At the 
same time, these versions give an acceptable description of the meson leptonic 
decay characteristics. 

PACS number(s): 11. 10. St, 12.39.Ki, 12.39.Pn 
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I. INTRODUCTION 



The Bethe-Salpeter (BS) equation provides a natural basis for the relativistic treatment 
of bound qq systems in the framework of the constituent quark model. However, due to 
the lack of the probability interpretation of the 4-dimensional (4D) BS amplitude as well as 
due to serious mathematical diseases which are inherent in the BS approach to the bound- 
state problem, various 3-dimensional (3D) reduction schemes of the original BS equation are 
usually used. As it is well known, the simplest version of this sort of reduction immediately 
arises if the kernel of the BS equation is taken in the instantaneous (static) approximation. 
As a result, the Salpeter equation is obtained. The Salpeter equation was used for the 
description of the bound qq system without further approximation in Refs. , whereas 

some additional approximations were made in Refs. [I3,n| . However, as it is well known, the 



Salpeter equation itself is not free from some drawbacks. Namely, it does not have a correct 
one-body limit (the Dirac equation) when the mass of one of the particles tends to infinity. 
From the general viewpoint, this property is expected to be important for the qq system 
with one heavy and one light quark. In order to avoid the above difficulty, in Refs. |15|JT6| 



the effective noninteracting 3D Green function for two fermions was chosen in the form that 
guarantees the correct one-body limit of 3D relativistic equations with the static BS kernel. 
These versions of the 3D equations will be referred hereafter as to the MW and CJ versions, 
respectively. A new version of the effective free propagator for two scalar particles which 
also possesses this property was suggested in Ref. [PT7|] . The effective 3D Green function for 
two noninteracting fermions can be constructed from this propagator in a standard way. 

Taking into account the fact that the relativistic effects are important for qq systems with 
quarks from light and light-heavy sectors, it seems interesting to carry out the investigation 
of this sort of systems in the framework of the above-mentioned different versions of 3D 
relativistic equations. This will allow one to shed light on the problem of ambiguity coming 
from the choice of a particular 3D reduction scheme of the BS equation, and to find the 
characteristics of the bound qq systems, which are more sensitive to this choice. For the 



meson mass spectrum this problem was addressed to in Refs. |T8|JT9|] where the MW and 
CJ versions of 3D relativistic equations together with the Salpeter equation (Sal. version) 
were considered in the configuration space to (partially) avoid the difficulties related to a 
highly singular behavior of the linear confining potential in the momentum space at the 
zero momentum transfer. The version of the 3D reduction of the BS equation suggested in 



Ref. [17] significantly differs from the MW and CJ versions and can be written down only 



in the momentum space. Consequently, it seems interesting to study together all versions in 
the momentum space and to investigate a wider class of characteristics of the bound systems, 
including the decay characteristics of the mesons, which are sensitive to the behavior of the 
meson wave functions, and the meson mass spectrum. These problems will form the subject 
of the present paper. 

The layout of the present paper is as follows. In Section II, we present different versions 
of the 3D bound-state equations in the momentum space, and perform the partial-wave 
expansion of the obtained equations. The numerical solution of these equations with the 
oscillator-type potential is considered in Section III. The general structure of the meson 
mass spectra obtained from the solution of these equations is discussed in detail. In Section 
IV, the leptonic decay characteristics of the pseudoscalar and vector mesons are calculated 
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using the wave functions obtained from the solution of these equations. 



II. THE RELATIVISTIC 3D EQUATIONS 

The relativistic 3D equations for the wave function of the bound qq systems, correspon- 
ding to the instantaneous kernel of the BS equation, i.e. when K(P;p,p') — > K st (p,p'), for 
all versions considered below in the c.m.f. can be written in a common form 

$m(p) = G 0eff (M,p) J |^ [iK st (p,p') = V(p,p')] $ m (p') (1) 

where M is the mass of the bound system, and the equal-time wave function <& M {p) is 
related to the BS amplitude $p(p) as 

**W) = f^*P=(M,5)<P) (2) 
The effective 3D Green function of two noninteracting-quark system Goe// is defined as 

G 0eff (M,p) = J ^[G 0eff (M-,p) = g 0eff (M-,p)(p' 1 +m 1 )(p' 2 + m 2 )] (3) 

Here g 0e ff(M;p) is the effective propagator of two scalar particles. The operator G 0e // is 
given in the form 

Goe ff (M,p) = £ £ D{ Z(M M ^ V) A » ^ 7?72 °' P S 1*1 (4) 

a 1= ±a 2 =± a \ lvl iP) 

where the projection operators A^ 1 " 2 '' are defined by 



2ui 

hiPi) = liiliPi) + ma®, Wi = (mf + ft 2 ) 1/2 (5) 
and the functions Z)( 0l02 )(M; p) and d(M;p) are given by the expressions (see Ref. [EH]) 

£)(aia 2 ) _ v / d = 1 

lu 1 + lu 2 - {oi\Ei + a 2 -E , 2) ' 

2 2 

+ £ 2 = M, E l -E 2 = mi M m2 = b o (MW version) (6) 

d = 2(wi + u^a, a = E?-Ui = [M 2 + 6g-2(a;i +w|)]/4 (CJ version) (7) 

Note that for the case of CJ version Eq. (|7|) is obtained from Eqs. (|3|) and (|j) by using 
the expression for g 0e ff(M]p) defined from the dispersion relation which guarantees the 
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elastic unitarity. The same condition is satisfied by the expression of go e ff(M;p) suggested 
in Ref. 0], (see formula (10) from this paper). According to this condition, the particles in 
the intermediate states are allowed to go off shell proportionally to their mass, so that when 
one of the particles becomes infinitely massive, it automatically keeps that particle fully 
on-mass-shell and the equation is reduced to the one-body equation. Using this expression 
for g Qe ff(M;p) in Eq. (|), we derive the expression for Goe// having the form (f|) where 



2y 



R — b , . 
+ [Ei + aiUi) - (E 2 + a 2 oj 2 ) 



d = 2RB, R=(b 2 -4y 2 a) 1/2 , B = 

2y 



■2y 
R-b 



■2y 



+ b 



+ a, 



b = M + b y, y= ^l^l ( 8 ) 
mi + m 2 



Hereafter this version will be referred to as the MNK version (our comments on Ref. [17 
see in Ref. pOf 



Using the properties of the projection operators A^ 1 " 2 -* and Eqs. (§)-©, the following 
system of equations can be derived from Eq. (|T]) 



[ M - + a 2 u 2 )]^r 2 \v) = 

= A^\M ]P )A^ ) ip-p)J^^V{p i p') £ £ ^ A \p') (9) 

a[=± a' 2 =± 

where &m \P) = A^ 1 (p, —p) &m(p) and the functions A^ aia2 \M] p) are given by 

A (±±) = ±1 a (±t) = M ( MW vers i on ) ( i0 ) 

UJi + UJ 2 

A(ma 2 ) M + {a x u x + a 2 u 2 ) 

A 1 1 2; = (CJ version) (11) 

2(a>i + uj 2 ) 

(MNK version) (12) 



-[M - (ai^i + a 2 u 2 )]— - — 



R — b . , . 

— h (Ei + a\U\) - (E 2 + a 2 u) 2 ) 



As to the Salpeter equation, it can be obtained from the MW version by putting = 

and l>£ T) = 0. 

Note that for quarks with the equal masses (mi = m 2 = m) and u = (m 2 + p 2 ) 1 ^ 2 , from 
Eqs. (0) and (0) we have 

A (±±) = ±1; A (±T) = M_ ( MW version) (13) 



2c; 
M 

Auj Auj 



A(±±) = M + 2u a (±t)=M (CJ version) (14) 
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One observes from Eqs. (p~3|) and ( |14j ) that in the limit of the equal-mass quarks the bound- 
state mass M enters multiplicatively in the coefficients in front of the mixed-energy com- 
ponents &m T \p) both in the l.h.s. and r.h.s. of Eq. (|9|). Consequently, dividing both 
sides of the equations for these components by M, one arrives at the (nondynamical) con- 
straints on all components of the wave function which must be considered together with the 
remaining two dynamical equations for the components ^m^Cp)- These equations for the 
bound state mass M are linear in the MW version and are nonlinear in the CJ version due 
to the fact that the r.h.s. of the equations depends on the value of M one is looking for. 
In the MNK version with an account of the property of the function R given by Eq. (||) 
lim R/y = lim M/y, from Eq. (O) we have 

mi-tm2 mi^m2 1 * 

A(±±) = M + 2u A (±T) = I (MNK version) (15) 

2,uj 2 

In this case, from Eq. (|9|) it follows that characteristic features of the bound-state 
equations remain unchanged - again it is the system of nonlinear equations in M and it 
includes all components of the wave function ^^ ia2 \p), as it is in the case for the quarks 
with nonequal masses. 

~ (a a ) 

Further, we write the unknown function "J^ 1 " 2 in Eq. ([|) in the form analogous to that 
used in Ref. ||, where the bound qq systems were studied in the framework of the Salpeter 
equation 

.(aia 2 )/->x _ Ar(aiaa)/ s ( \ „ / ' \ , >l<* 2 ). 



•ff-'OO = jvjf-o,) ( ai(<fi?r)/( ; + ^ j s { _ a2{9 ^ 2 + a2m2) ) 

(16) 



where 



N- 



(aiQ2) _ fun + gimi \ V 2 (u 2 + « 2 m 2 \ V 2 _ A A ai ) ( , AT (a 2 ) 
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( P ) = { 1 2 j; i ) [ z 2 ;; z ) ^Nrmnp) (17) 



Then, if the qq interaction potential operator V(p,p') is taken in the form || 

V&P,P') = lhlV og {p-p') + [x 7 ?7 2 ° + (1 - x)hI 2 ]V c (p-p'), (0 < x < 1), (18) 
the following system of equations for the Pauli 2 <g> 2 wave functions ^m" 2 ' can be derived 



[M - (a^ + a 2 u 2 )} X { M ia2 \p) 

d 3 p 

W) 



where 



)5 1 (p,p ,ai,a 2 ) 



+ ^(x;p-p )5 2 1 2 (p,p ,a u a 2 ) 



N[T 2 \p>) (20) 
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j^(a 1 a 2 d 1 a 2 ) _ ^ _|_ 



a.\a.id x cl 2 (°lP ) (^2P ) (glP Q (ggP Q 

(c^i + Qimi)(w 2 + a 2 ^2)(^i + «'i^i)(^2 + ®' 2 m 2 ) 



(21) 



^(aia 2 o 1 a 2 ) 



a 1 a' 1 (o ! 1 p)(o ! 1 p ') a2a' 2 (a 2 p)(a 2 p ' ') 



(ooi + a- l rn- i ){u l + a^mi) (tu 2 + "2^2) (^2 + «2 m 2) 



^(x;p-p ') = K*(p-p') + (2x- l)V c (p-p ') 
Now using the partial-wave expansion 

X$T°GP)= E (P I LSJM^R^ip), (P = l) 

LSJMj P/ 



(22) 
(23) 

(24) 



where 



V(p-p') = (2vr) 3 E V L { P ,p')(p\LSJMj)(LSJMj\p) 

LSJMj 



V (p,p) = - / j L (pr)j L (pr)r 2 dr 

7T JO 



(25) 



(26) 



with j'l(x) being the spherical Bessel function, the following system of equations can be 
obtained from Eq. (|19D for the radial wave functions R^sj^ip) 

[M - {a lUl + a 2 u 2 )\R^?{p) = A^(M;p)N^ Q2 \p) x 



x 



roc , , 

E / p 2d p 

Jo 



i + a&'S' 4a4) (p,i/)V J (i;p,p) + 



+ a^ Q2 ' QlQ2) (p,p>i" J ) (x;p,p) 



(0102,0102)/ '\T/ / '\ 

a e (p,p )l/ eJ (x;p,p ) 



x 



[ M - (a^! + a 2 w 2 )]fl§±nJ(p) = A^{M;p)N& a2) (p) x 

r r JO L . 



( Q l a 2) 
JillJ 



(?') + 



+ 



(ai02,Oi0 2 ) 
l 12(g>12 



2J+ 1 



Vej(i;p,p') 



(27) 



where 
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e 

4*^W) = a?(p)a?(p) , a?(p) = ^ P (28) 



and 



V, 



$ & p, P') = [ ( J + x ) 7 J - 1 (x; p, p') + ( J + 1 ) V J+1 (x; p, p')] 

V eJ (x;p,p') = ^ J ^j^\ v J - l {x-p,p') - V J+1 (x;p,p')} (29) 

Vj ±11 j{1;p,p') = * [V J±1 (l;p,p')+iJ{J + l)V J ^{l;p,p')) 

V(x;p,p') = V og (x;p,p') + (2x - l)V c (x;p,p') 

The main purpose of the present study is to carry out the comparative qualitative ana- 
lysis of the different versions of 3D relativistic equations (p7|), addressing the question of 
existence of stable solutions for different values of the scalar-vector mixing parameter x in 
the confining part of the potential (Eq. (|T8"D). Also, we investigate the general structure 
of the meson mass spectrum and calculate the leptonic decay characteristics, namely, the 
pseudoscalar decay constant fp(P — > \w) and the vector meson decay width T(V — ► e~e + ). 
For these reasons, at the first stage we neglect in ([Tj|) the one-gluon exchange potential. 
A full analysis of the problem will be made in further publications. Further, according to 
Ref. PH , we use the oscillator form for the confining potential V c (r) which is a simplified, 
but justified form of a more general potential used in Ref. || (at least for the quarks from 
the light and light-heavy sectors, which are considered in the present paper). Namely, we 
take 

K(r) = ^K 2 )(^r 2 -V ) (30) 

^im 2 2 12-k ( Q 2 \- x 

A*i2 = , m 12 = m 1 + m 2 , a s {Q ) = — — - — hi-^ 

m 12 33 — 2rif \ A z J 

where n/ denotes the number of flavors. This potential in the momentum space corresponds 
to the operator (see (p6|)) 



K L (P ? P) = ~®s{m 2 12 ) 



VuuL d 2 2_^_ _ L(L + 1) 
2 dp z p dp p z 



s(p-p) (31) 



ppi 



Now the system of equations (P7|) can be reduced to the system of equations with the 
following structure: 



[ M - (atUJx + a 2 u 2 )]R^\p) = -^a s (m 2 12 )A^\M;p) £ I V 



B { « ia2 ' a ^\p) + 
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£>B (P) + ( 2x ~ l ) D A (P) ~ 



) (aia 2 ,a 1 a 2 ) . 



p2 ' 



-(2x 



'(p)J(J + l) + (2x-l)4 



j(j + r 



p- 



[M - (a lWl + a 2 c, 2 )]4^ 2 j(p) = -|a s (m? 2 )^ Q2 )(M;p) £ 



5, 



(aia2,a 1 a 2 ) 



(p) + 



+(2x- 1)4 



(aia 2 ,a 1 a 2 ) 



(p) 



1? 



(aia 2 ) 
J±11J 



(P) + 



J- / -,(aia2,a 1 a 2 ) / « / - - . ! 



— r B 



(p) J(J + 1)+1± 



± 



2 J + 1 



+4 T J i J t! ; f 4 aia2,a ' ia2) (p) - 4 Q1Q2 '° lQ2) (p) XTx^p) 



4J( J + 1) R (ttia 2 ,a' 1 4), 



r v ~ ' _, v ' 2j + i, 

+ (2* - l)A { ^ a ' A \p)J(J + l))K±n](p) + 



'(P) + 



p 2 (2J+1) 



(32) 



where A^^'^^anciBi" 102 '" 1 " 2 '' are a given functions of p, and Z)^^ 2 '" 1 " 2 "' arc certain 



e e 
second order differential operators in p. 



III. MESON MASS SPECTRUM 



In order to solve the system of equations (B^) for the bound state mass M, the unknown 
radial wave functions R^sjip) are expanded in the basis of the radial wave functions 
Rul(p) being the solutions of the Shrodinger radial equation with the oscillator potential in 
the momentum space Q3"T|), as it was done in Refs. 



R?S?\P) = (2M(2n)^ 2 CfsZ ] RnL{p) 

?1=0 



(33) 



where the multiplier (2M(27r) 3 ) 1/ ' 2 is introduced for appropriate normalization of the wave 
function and 



Rul(p) =P 
P 

z = — , 
Po 



2r (n + i + |) V /'?^pK) iFi( _ K!L + i;z2) ^ iW 



r(n + i) 



r(L + |) 

l/2\ 1/2 



(34) 



Substituting the expression (|33|) into the system of differential equations (0), the following 
system of algebraic equations for the coefficients C^sfn can be obtained 



S 



Me*-' = e e (35) 



It is necessary to note here that the matrix H a/3 explicitly depends (except the Sal. 
version) on the meson mass M we are looking for. Consequently, the system of equations 
(35) is nonlinear in M. Note also that for the quarks with equal masses (mi = m 2 = m) part 
of the equations from the system (|35|), corresponding to (ai« 2 ) = (±T) f° r the MW and 
CJ versions, transforms into the nondynamical constraints between all coefficients Clsj^ 
which should be considered together with the remaining dynamical equations corresponding 

to (tti« 2 ) = (±±)- 

The numerical algorithm for the solution of the system of nonlinear equations (|35D 
in the case of nonequal mass quarks was discussed in Ref. [^(J where the systems us 
(^i, 1 Pi, 3 P , 3 Pi, 3 P2, 1 P2, 3 Pi, 3 P 3 ), cu and cs (^o, 3 -V Pi, 3 P 2 ) were considered. In brief, 
the infinite set of equations (^) is truncated at some fixed value n = N max and the eigen- 
value problem is solved for the finite-range matrix H. Increasing then N max , one checks the 
stability of the resulting spectrum with respect to the variation of N max . Since the r.h.s. of 
Eq. ([HJ) depends on M, the solutions are obtained iteratively, starting from some value of 
M. In Ref. pOfl , the existence of stable solutions of the system of equations (^) was inves- 
tigated for different values of the mixing parameter x in the confining potential flTBp. It was 
found that the existence/nonexistence of stable solutions of Eq. (^) critically depended on 
the value of x, on the value of confining interaction strength parameter t^o (|3"TD, and on the 
particular state under consideration. This dependence is different for the different versions of 
3D reduction of the BS equation. The instability is primarily caused by the presence of the 
mixed (H — , — h) energy components of the wave function in the equations for the qq bound 
system. Namely, for the parameter u =710 MeV that leads to a reasonable description of 
the meson mass spectrum in the framework of the Salpeter equation ||, stable solutions 
for CJ, MNK and Sal. versions simultaneously exist for the values of the parameter x from 
the interval 0.3 < x < 0.6 — 0.7. For the MW version, in order to provide the existence 
of stable solutions in the same interval of x, u>o must be set to a smaller value (450 MeV). 
However, in this case the values of masses for all states under consideration turn out to be 
smaller than the experimental ones. Keeping in mind that the calculated values of masses 
will further decrease after adding the one-gluon interaction potential, we conclude that the 
MW version seems to poorly describe the meson mass spectrum. For this reason, along with 
the Sal. version, below we consider only the CJ and MNK versions, both having a similar 
theoretical foundation: the effective Green function (^) in these versions is constructed from 
the elastic unitarity condition. 

The results of calculations are given in Figs. 1,2,3, from which one can see that the 
level ordering is similar for all three versions under consideration. Further, at x = 0.5 the 
states 3 Po, 3 Pi, 1 Pi are degenerate and spin-orbit splitting exists only between the degenerate 
3 Po, 3 Pi states and the 3 P 2 state. For x ^ 0.5 this degeneration is removed and the calculated 
level ordering agrees with the experiment for the value x = 0.3, except the 3 P 2 state. For 
the D - states (us), experimentally there is degeneration between 1 P 2 , 3 P 3 states. In our 
calculations we do not have this degeneration, but for the MNK and CJ versions at x = 0.3 
the splitting is very small and increases for x = 0.5 and x = 0.7. Note, however, that only 
for these values of x the sequence of the 3 Pi state and other two D-states agrees with the 
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experiment. For the qq states with the quarks from light-heavy sectors (cu, cs) the spin- 
dependence of the energy levels for all values of x is much weaker than the experimental 
one, but at the same time for x > 0.5 the level ordering agrees with experimental data. 

As to the pseudoscalar qq systems (the 1 Sq state) with the quarks from the light sector 
(us), the calculated masses in the model under consideration are much larger than the 
experimental ones, as demonstrated in Fig 1. This might serve as an indication of the fact 
that if the constituent quark model is used for the description of this sort of systems, the 
chiral symmetry breaking effect should be taken into account, at least in a phenomenological 
manner, e.g. by the inclusion of the t'Hooft interaction in the kernel of the Salpeter equation 
(see Ref. ||). In order to take into account a full content of global QCD symmetries in a 
systematic way, a coupled set of Dyson-Schwinger and BS equations should be considered 
(see, e.g. 0). 

Note that the number of terms (N max ) in the expansion (^), which is necessary to get 
stable solutions of the system of equations fl3"5]), varies with the constituent quark masses, 
with the value of the mixing parameter x, with the state under consideration and is differ- 
ent for the CJ, MNK and Sal. versions. Namely, when the quark masses increase, N max 
decreases. The convergence of the numerical procedure used in the calculations is better for 
x < 0.5 and worse for x > 0.5. For all values of x the convergence is better for the Sal. 
version than for the CJ and MNK versions. 

We have also calculated the mass spectrum of qq systems with the equal quark masses 
from the light quark sector (uu, ss). This problem was the subject of our primary interest in 
the present study. The results of calculations are shown in Figs. 4 and 5. Note that for this 
sort of systems the convergence of numerical algorithm used in the calculations appears to 
be not so good for the values of the parameters u =710 MeV and x > 0.5. The convergence 
becomes better for smaller values of lvq. Namely, for the (uu) system for c<Jo=710 MeV and 
x = 0.6, stable solutions in the MNK version for the 3 S\ state do not exist. In the CJ 
version such a situation holds for other states ( 3 P2, 3 -Di, 3 -D3) as well. For smaller values of 
the potential strength parameter, e.g. u;o=550 MeV, the stable solutions exist for all states 
(just these results are shown in fig. 4). Further, in this case the sequence of the energy 
levels corresponding to the 3 Pj states (the spin-orbit splitting) agrees with the experiment 
at x > 0.5 in all versions, and in the 3 Di and 3 D 3 states the agreement appears to occur at 
x < 0.5. 

Consequently, on the basis of the above analysis one can conclude that none of the 3D 
equations with the simple oscillator kernel considered in the present paper, does simultane- 
ously describe even general features of the mass spectrum of all qq systems under study. 



IV. SOME DECAY CHARACTERISTICS OF MESONS 



For the investigation of the meson decay properties, the normalization condition for the 
wave function $m(p) — J2 ®m (p) is needed. For the Salpeter wave function this 

(0102) 

condition is well known |23| 



d 3 p 
J2n)~ 3 



2M 



(36) 
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As to the MW, CJ and MNK versions, the normalization condition for the corresponding 
wave functions can be obtained with the use of the fact that the effective Green operators 
in Eq. can be inverted. As a result, the corresponding full 3D Green operators Go e ff 
can be defined analogously to the 4D case 



G:} f (M;p,p') = (27r) 3 5^(p-p')G 0e 1 ff (M;p) - V(p,p ') 

Since V{p,p') does not depend on M, the normalization condition reads 
d 3 p 



(2vr) : 



3> m (p) = $m(p) 7i U 7 2 U 



DM 



$ M (p) = 2M. 



(37) 



(38) 



Using Eqs. 



from Eq. ( p8|) one obtains 



E 



~<\ aia2 \p)flr 2 \M-,p)&M ia2 \p) = 2M 



where 



(2tt) 3 

( ai a 2 ) _ ogEx + a 2 E 2 



< 12 



M 



Aaxa 2 ) _ u l + u 2 a\U) 1 E 2 + a 2 UJ 2 Ei 
J 12 ~ 



M (Ei + ociUi) (E 2 + a 2 oj 2 ) 



f (0102) 
J 12 



M 



5(l-y 2 ) + 



M 2 



- 2 



(MW version) 

(CJ version) 
R-M 



(39) 

(40) 
(41) 
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+ (aicoi — a 2 cu 2 



M + 



a\iO\ + a2W2 



M 2 



- 2 



2y 



R-M\- 



+ 



i?-M M 



2 V 2y 

(MNK version) 



(42) 



Note that for mi = m 2 = m the normalization condition (^) for the MW version is reduced 
to fl3~6|) which can be written in the form of Eq. ( |39"D where 



f (aia 2 ) _ a l + Q 2 
Jl2 — o 



(Sal. version) 



(43) 



From the normalization condition ( p9| ) for the wave function given by Eqs. (|T6|), flTTD one 
obtains the normalization condition for the wave functions R^sj^Ip) with the use of the 
partial-wave expansion (f24|) 



o (2vr) 3 ^ 2 



R 



(ai«2) 
LSJ 



(P) 



n 2 



2M 



(44) 



The functions /{ 2 ^ (40,41) have second order poles at p = p s , where 

a{p s ) = Ef - a$(p.) = 0, p s = i(M 2 + b 2 - 2{m\ + m 2 )) 1/2 



(45) 



11 



The functions /{J turn out to be finite (f\2 (p s ) is apparently finite). Consequently, 
the normalization condition (pH ) for the CJ and MNK versions involves a singular integral 
of the type 



I(x ) 



f(x)dx 



(46) 



/o (x — x ) 2 

which taking account of the conditions /(0) = = /(oo) valid in our case, can be regularized 



as 



f(x)dx = f 2x ° [fjx) - f'(xo)]dx f 00 f(x)dx 

(x — Xq) 2 JO (x — Xq) J2x X — Xq 



(47) 



Now we can calculate the pseudoscalar (S = L = J = 0) decay constant fp(P — > \iv) 
and the leptonic decay width of the vector (/ = 0, S = J = 1) meson T(V — > e~e + ) 
(the corresponding decay constant is denoted by /y). In these calculations, instead of the 
functions ^^ lCl2 \p) ( pTC|) given as a column with the components $^" lQ!2 \ $„" ia2 , $^ ia2 \ 

bb 5 

Ref. IE3I 



$^ 1 " 2; , it is convenient to introduce the wave function \I/( aia2 ) written in the form (see 



,7,(0:10:2) _ ^no ^ab //^ _ ,-^.2 0\ 



$r 1U2J $L 

ba bb 



(48) 



where C is the charge conjugation operator. Then, the decay constants fp and fy are given 
by the expressions 



SuaMfp = V3Tr 



*ooo(r = 0)7^(1-75) 



(49) 



f v {\) = V3Tr 



*oua(^= 0)7 a 



e"(A = 0,±1) 



(50) 



Here, the factor \/3 stems from the color part of the wave functions, £ M (A) is the polarization 
vector of the meson and 



LSJMj 



(f = 0) 



<i 3 p 



LSJM 



,W = E * 



(aia 2 ) 
LSJM j \P ) 



Using Eqs. (0), (0), (0), (g§) and (0), from Eqs. (|p]) and (|5Jj) one obtains 



/i 



M 7o (2tt) 3 ^ 



0'2 



E N^\p)Nt\p)-a 1 a 2 Nt ai) (p)Nt a2 \p) 



R ( Z° 2 \P) (52) 



(51) 



MA) 



-V24 



p 2 dp 



Ek w (p)ivf 2, (p) + 



+ 



o (2tt)3^ 

^^ ( - Qi) (p)iv 2 ( - a2) (p)]4?r 2) (p)}5 



(53) 
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For a given fy ([53]) the leptonic decay width of the vector meson is given by 



<ff 1 
4ttM 3 3 



E IMA) 



A=0,± 



where 



a 



1/111 



e// 137V2'18'9 



for q°,uj and mesons, respectively. 



Finally, using Eqs. fl3|), (£§, (§2]), (§3|) and fl5|), we obtain 



/i 



E 



7TVM Q, iaa ^0 



« 1 a 2 iV 1 ( - ai) (p ^)iVr Q2; (po^) 



z dz 



Ni ai) (p ,z)N^\p ,z)- 



(-£* 2 )/ 



-"-000 l z J 



(54) 



(55) 



iyr^o^)ivr 2j (po^) + 

+ ^iv(--)(p 0) ,)iv(--) ( p 0) ,) 

where the functions 

oo 

R { ls?\z) = E ^j(«i« 2 ) j R.l(^) (57) 

n=0 

satisfy the normalization condition 

2 = 1 (58) 

The results of numerical calculations of the quantities fp and T defined by Eqs. ( |55| ) and 
(0), are given in Table I. We see from Table I that the calculated values of fp in the MNK 
and CJ versions, as a rule, are smaller than in the Sal. version and this fact is related to the 
presence of the contributions from the "H — " and " — h" components of the wave function 

(contribution from the " " component is negligibly small). Further, the calculated value 

of fp is larger in the CJ version than in the MNK version. The calculated value of the 
quantity T(V — > e~e + ) weakly depends on the particular choice of the 3D reduction scheme 
of the BS equation. With the increase of the mixing parameter x both the quantities fp and 
T(V — > e~e + ) slightly increase. The calculated values of fp and T{V — ■> e~e + ) are smaller 
than the experimental ones. 



ViV ^ e 6 )= (27W' ?Jo ZdZ 



tCr\z) | 2 (56) 



E 



z 2 dzft a2) (M;p , 



R { ls?\z) 



On the basis of the analysis of the different versions of the 3D reductions of the bound 
state BS equation carried out in the present paper, one arrives at the following conclusions: 

The existence/no nexistence of stable solutions of the 3D bound-state equations criti- 
cally depends on the value of the scalar- vector mixing parameter x. For all 3D versions 
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(Sal, MNK, CJ) stable solutions coexist for the value of x from a rather restricted interval 
0.3 < x < 0.6 — 0.7. The level ordering in the mass spectrum is similar for all versions 
under consideration. However, the sequence of the calculated energy levels agrees with the 
experiment for some states at x < 0.5 and for other states at x > 0.5. Consequently, a si- 
multaneous description of even general features of the meson whole mass spectrum turns out 
not to be possible for a given value of x from the above-mentioned interval. It is interesting 
to investigate the dependence of this results on the form of the confining potential. Also, 
it is interesting to study how it changes when the one-gluon exchange potential is taken 
into account. This aspect of the problem will be considered at the next step of our inves- 
tigation. Further, we plan to include the 't Hooft effective interaction in our potential in 
order to (phenomenologically) take into account the effect of spontaneous breaking of chiral 
symmetry which is important in the pseudoscalar ( l So) sector of the constituent model. 

The calculated leptonic decay characteristics of mesons are quite insensitive to the par- 
ticular 3D reduction scheme used and give an acceptable description of experimental data. 
In future, we also plan to study the validity of this conclusion for a wider class of realistic 
interquark potentials. 
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TABLE CAPTIONS 



Table I 

The pseudoscalar decay constant fp(P — > iiv) (in MeV) and the leptonic decay width 
T(V — > e~e + ) (in KeV) with the allowance for only (++) and all components of the wave 
function. 



FIGURE CAPTIONS 

Fig. 1 

The mass spectrum (in GeV) of the us system for the different 3D equations and different 
values of the scalar-vector mixing parameter x. The multiplicity of degenerate levels is 
indicated by the number inside the dash (u = 710 MeV). 

Fig. 2 

The same as in Fig. 1 for the cu system. 
Fig. 3 

The same as in Fig. 1 for the cs system. 
Fig. 4 

The mass spectrum (in GeV) of the uu system for the different 3D equations and different 
values of the scalar-vector mixing parameter x. The multiplicity of degenerate levels is 
indicated by the number inside the dash. The parameter lu for the 3 Si state in the MNK 
version and for all states in the CJ version at x — 0.6 is fixed at the value 550 MeV, 
otherwise oo = 710 MeV. 

Fig. 5 

The same as in Fig. 1 for the ss system. 
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Table I 



Decay characteristics 




fp(P- 


-> nv) 






r(v -» 


e-e+) 




Meson 


D(cd) 


n ( 


cs) 


g(uu) 


ip(ss) 


Versions 


0:1^2 


x=0.3 


£=0.5 


x — u.o 


X — u.o 


x=0.3 


x=0.5 


x=0.3 


a;=0.5 


Sal 


++ 


149. 


156. 


177. 


183. 


3.94 


4.31 


0.84 


0.90 




all 


148. 


155. 


176. 


182. 


3.96 


4.38 


0.84 


0.90 


MNK 


++ 


148. 


155. 


176. 


181. 


4.44 


4.85 


0.89 


0.95 


» 


all 


141. 


145. 


168. 


172. 


4.11 


4.20 


0.85 


0.89 


CJ 


++ 


150. 


158. 


179. 


186. 


3.98 


4.45 


0.86 


0.94 




all 


149. 


152. 


176. 


180. 


3.91 


4.23 


0.85 


0.90 


Expt. 


< 220 


170 -: 


- 180 


6.8 z 


t0.3 


1.37 =1 


= 0.05 
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